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Abstract 

In this short reappraisal of spontaneous lepton number violation in a supersymmetric scenario 
, implemented through singlet sneutrino vacuum expectation value (VEV), we contribute with two 

t-H ■ new things in the context where the lepton number symmetry is global: (i) provide explicit ex- 

pressions of i?-parity violating couplings in terms of the neutrino Yukawa couplings and the singlet 
sneutrino VEV, and (ii) estimate the limit on this VEV using the current knowledge of the light 
neutrino mass and the astrophysical constraint on the Majoron-electron coupling. Besides, we put 
updated constraints on the VEV and Yukawa coupling of the singlet superfield when lepton number 
| is gauged. 
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■ 1 Introduction 

O It is well-known that if the supersymmetric partners of all standard model (SM) particles are intro- 
duced in a theory and one constructs the most general Lagrangian that is invariant under supersym- 
metry and the SM gauge symmetry, the Lagrangian contains terms which violate both lepton number 
(£) and baryon number {B). The explicit C- and B- violating parts appear in the superpotential as: 

KbcL a L b E e + X ab c L aQbD c + X'^ bc U a D b D c + ^ ^ a L a H u . (1) 

■ abc abc abc a 

Above, all superfields are left-chiral, and the subscripts a, b, c are generation indices on lepton doublet 
fields L, quark doublet fields Q, and SU(2)-singlet charged fields E, U and D, in obvious notations. 
The hat on a superfield means that the left-chiral fermion part of that superfield is the antiparticle 
of fermion whose name is suggested in the letter denoting the superfield. These terms also violate 
imparity, defined by (— i) 3B +£+ 2i5 5 where S is the spin of the particle. The antisymmetry in the first 
two generation indices of A and in the last two indices of A" suggests \ abc = — A& ac and A" 6c = — A" cb . 
Clearly, there are 9 A-type, 27 A'-type and 9 A"-type trilinear, plus 3 //-type bilinear, i?-parity violating 
(RPV) couplings. These 48 new couplings add further twists and complications to phenomenology 
and bring in more uncertainty to theoretical predictions. 

Our aim in this paper is to explore a restrictive scenario in which there will be much fewer RPV 
couplings, thus offering more predictivity. For definiteness, we assume that these couplings are gener- 
ated by spontaneous C violation. This immediately rules out the ,6-violating A"-type couplings. Now, 
we recall that Aulakh and Mohapatra pQ were the first to have implemented the idea of spontaneous 
violation of C = 1 global lepton number in a supersymmetric context through the VEV of the sneu- 
trino component of a lepton doublet superfield. A testable feature of this realisation was a photino 
mediated contribution to neutrinoless double beta decay. Neutrinos were predicted to be massless at 
tree level, with a suggestion that very small masses (m v ~ (10 -5 — 10 -8 ) eV) are induced at one- loop 
order through a combination of supersymmetry breaking and lepton number violation. A follow-up 
study [2] revealed that if supersymmetry breaking terms include trilinear scalar couplings and gaugino 
Majorana masses, neutrino mass would be generated at the tree level itself, with a special property that 
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even with three generations only one non-vanishing mass eigenvalue would emerge at tree level. The 
mechanism of pQ would indeed induce two other small masses at one- loop order. Other implications 
of this scenario were studied in the context of matter-enhanced solar neutrino oscillation [3HH] . 

An important feature of spontaneous C violation is the existence of Majoron (J), which is a physical 
massless Nambu-Goldstone boson arising from the imaginary part of the sneutrino field that acquires 
a VEV. The mass of the real scalar (/?), associated with the Majoron, in the doublet Majoron scenario 
turns out to be very small leading to unacceptably large Z — > pj decay, which is ruled out by the LEP 
data on Z boson decay width. In fact, gauge non-singlet Majoron models are all strongly disfavored 
by electroweak precision measurements OH]- Subsequently, singlet Majoron scenarios were proposed 
in the supersymmetric context. In some of these models, lepton number was spontaneously broken by 
the VEV of a field carrying one unit of lepton number [HHlOj , and in some others by a field carrying two 
units of the lepton number [11], like the non-supersymmetric models of spontaneous lepton number 
violation |12[I13], Since the first kind, i.e. AC = 1 violation, is a speciality of supersymmetric models 
that non-supersymmetric models do not have, we take it up for our work here. 

In this context, we present explicit expressions of X a bc, X' abc anc ^ A*a couplings in terms of the 
Yukawa couplings of the general superpotential and the singlet neutrino VEV. Furthermore, we provide 
new bounds on the sneutrino VEV and the generic neutrino Yukawa couplings from astrophysical 
considerations of stellar energy loss. 

2 imparity violating couplings 

Our model has, apart from the superfields in the MSSM, gauge singlet superfield N a , one for each 
generation. These fields carry lepton number C = — 1, and hence their VEVs (V a ) would spontaneously 
break lepton number. The superpotential of this model can be written as 

^ = z ^h i ^Q a U b H u + Y,h i $QaD b H d + Y,h a 2L a E b H d + ^ (2) 

ab ab ab ab 

We assume that the VEVs of N a are generated by the same mechanism as in [8l-H0j which require 
the existence of two more gauge singlet superfields (with C = +1 and 0, respectively). Though we 
implicitly acknowledge their existence we do not explicitly display how those two additional singlets 
appear in the superpotential, whose raison d'etre is to provide the V^'s. Beside that, Eq. (|2J) is 
the most general gauge invariant superpotential that also conserves lepton number before the scalar 
components of N a go to the vacuum. It is not difficult to realize that the VEVs V a by themselves cannot 
break super symmetry. Also note that apart from the term containing the N fields, the superpotential 
corresponds exactly to that of i?-parity conserving minimal super symmetry. 

The /^-violating operators are generated as soon as the VEVs V a are induced. Fig. [Tk. will generate 
the A terms, whereas Fig. [lb will generate A' terms of Eq. ([I]). The important point is that, these 
couplings will now be determined by Yukawa couplings and the sneutrino VEVs. It should be noted 
that when N acquires a VEV, the internal line in Fig. [T] is necessarily Higgsino. If we assume that 
the Higgsino mass parameter p is a few hundred GeV or more, one can effectively write a contact 
interaction from Fig. [TJ When the singlet sneutrinos acquire VEVs, we obtain 

- E! (3) 

d M 

d v 
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Figure 1: Generation of LLE and LQD operators when the scalar part of N acquires a VEV. The thick lines 
denote superfields. Generation indices have been suppressed. The blob in the middle of the internal line implies 
that it involves the fi-term of the superpotential. 

The two terms in case of A a f, c arise depending on whether L a or Lf, appears in the same vertex with 
the Nd superfield in Fig. [Th. The bilinear lepton number violating terms, shown in Eq. (pTJ) , also arise 
in this model from the terms of the superpotential in Eq. ([2|) when the scalar component of N 
acquires a VEV: 

*. = £*2°^ (5) 

b 

The origin of the relative minus sign between the two terms in Eq. (j3J) can be understood by 
keeping the SU(2) indices. If we denote the SU(2) index carried by L a and by a and (3 respectively, 
and put the SU(2) indices 7 and 5 on the internal H u and Hd superfield lines in Fig. [TJ then the 
diagram with L a coupling directly to Nd will contain the SU(2) factor 

(6) 

whereas the other diagram, obtained by interchanging L a and L&, should contain 

EfiryEySEaS = £a/3 • (7) 

Hence the minus sign in Eq. (jSJ), which makes the coupling antisymmetric in the indices a, b. 

We can now estimate how many unknown parameters are present in the /^-violating sector. Without 
any loss of generality, the couplings and can be taken diagonal in the generation indices, and 
they can be made real. In this case, these couplings are proportional to the masses of the down- 
type quarks and charged leptons, and are therefore known. The couplings are irrelevant for our 
discussion since they do not appear in the expressions of Eqs. ([3]), (j4|) and These will contain 
the up-type quark masses and the parameters of the Cabibbo-Kobayashi-Maskawa matrix. For our 
purpose, the relevant unknown parameters appear from h^ N \ and they are 9 in number. Besides, 
there are the three VEVs V a . To be more precise, there is actually only one independent VEV of the 
singlet sneutrino fields, since we can always make suitable linear combinations of the three N fields 
leading to the occurrence of a single VEV. This makes a total count of 10, instead of the 39 /^-violating 
parameters appearing in Eq. ([T|). 

3 Sneutrino VEV and Majoron-electron coupling 

The Majoron-electron coupling arises both at tree level and at one loop order from different inter- 
actions, which could be of similar magnitude. The tree contribution originates from the fact that a 
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Figure 2: Diagrammatic way of seeing how a doublet sneutrino gets a VEV. The dashed lines denote the 
neutral scalar fields contained in the supermultiplets, and the cross-hatched blobs denote their VEVs. 

non-vanishing V, the generic VEV of a N field, must accompany a non-vanishing vl, the VEV of a 
doublet sneutrino. This can be seen most easily from the fact that in the scalar potential of the model, 
the F-term of the H u field contains a term of the form fih^ LNH d , where in this expression only the 
scalar components of each superfield is implied and generation indices are suppressed. The diagram 
in Fig. [2] now clearly shows that the magnitude of should be given by 

„ ^Vv d 

where h^ N ' is the generic Yukawa coupling involving the N fields, Vd is the Higgs VEV contained in Hd, 
and Ms is a generic doublet sneutrino mass. The above expression can also be appreciated directly 
from potential minimization. In the tadpole equation dV/dvi = 0, the trilinear term fih^LNHl will 
provide a contribution fih^ N ^Vvd for the left-hand side, while the soft mass term MgL^L will yield 
M^vl- The minimization condition thus gives Eq. (jSJ). 

The non-zero value of the doublet sneutrino VEV induces, from the supergraph shown in Fig. [1^, 
a tree-level electron-Majoron coupling. Assuming, for the sake of illustration, that the soft mass of 
Hd is of the same order as fj,, this coupling is given by 



£j*^(h iN) ) 2 Vm e , (9) 



where the factor m e , equal to times Vd, ensures a chirality-flipping coupling. The loop induced 
contribution to the electron-Majoron coupling arises from the diagrams shown in Fig. [3l when one of 
the external TV-lines obtains a VEV. A rough estimate of the coupling thus generated is of the order 




Figure 3: Effective operators that give rise to a coupling between the electron and the Majoron when one of 
the external scalar lines goes to the vacuum. 
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where Mq is the heaviest mass in the diagram, either of the Z boson or of the neutralino (through its 
Higgsino component). The magnitude of the tree and loop contributions to g ee j could be of the same 
order, or one may dominate over the other, depending on the magnitude of the parameters Ms and 
Mo- A cancellation between the two contributions is unlikely and we avoid any such fine-tuning. 



4 Combined astrophysical and neutrino mass constraints 

There are stringent astrophysical constraints on any Majoron model. Majoron emission leads to stellar 
energy loss, and singlet Majorons may be emitted via Compton-like processes 7 + e — > e + J. The 
allowed leaking of stellar energy can be translated into a bound on the singlet sneutrino VEV. In fact, 
it turns out that Majoron coupling to the electron should be less than about 1CP 10 |14H17| from a 
study of the main sequence stars. Putting Ms ~ 100 GeV in the tree level contribution Eq. Q, the 
Majoron emission bound implies 

(/i (A 0)V<2MeV, (11) 

while from the loop contribution Eq. (|10p . for Mq = 100 GeV, the bound turns out to be 

(/i (A 0)V<lGeV. (12) 

Although the bounds in Eqs. (jlip and (|12p are quite different, both are independently significant, as 
the scalar mass (Ms) involved in Eq. ([9]) can be much larger than the neutralino (or, the Z boson) 
mass (Mo) in Eq. (|10p . The above bounds can be more stringent if, instead of main sequence stars, 
we use red giant stars, which give g ee j < 3 • 10 -13 [H]. However, we use the constraints from main 
sequence stars which are more reliably understood. 

It is interesting to note from Eqs. ([9]) and ()10p that the electron coupling to the Majoron is directly 
proportional to the lepton number breaking VEV. In contrast, the same coupling is inversely propor- 
tional to the lepton number breaking VEV in the non-super symmetric singlet Majoron model |12pi3j. 
The reason for the difference is that in the non-super symmetric case, where lepton number symmetry 
is broken by the VEV of a scalar field carrying two units of lepton number, heavy singlet neutrinos 
whose mass is of the same order as their VEVs float in the loop causing propagator suppression. In 
our case, V appears only in the numerator when an N is replaced by its VEV. 

We now discuss how neutrino mass is generated in our scenario. Fig. U] is a diagrammatic repre- 
sentation of AC = 2 Major ana mass generation at tree level. It gives 

m ab oc hVphUpVM. (13) 

Since the mass matrix is of rank one, only one non-vanishing eigenvalue will emerge. This is not 
surprising, as one can always perform a basis rotation in flavor space to put the VEV only along one 
direction. In fact, what we discussed is nothing but the neutrino mass generation through bilinear RPV 
couplings. Accurate expressions of the tree level neutrino mass induced by bilinear RPV couplings 
can be found, for example, in [Sj. For our purpose, it is enough to use the approximate expression of 
neutrino mass suggested by Fig. HI 

^~f 2 (" w ) 2 1?r' (14 > 
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Figure 4: Schematic diagram showing how a neutrino acquires a tree level mass in our model after the N 
fields and the neutral H u acquires VEVs. 

where v u denotes the VEV of the scalar component of H u , and My is a generic neutralino mass 
capturing the effects of the Zino and Higgsino propagators in Fig. HI where it is implicitly assumed that 
~ fi. Even though we do not yet precisely know the absolute magnitude of neutrino mass, we make 
a reasonable guess by putting m u = 0.1 eV in Eq. (fT4"|) , We further assume that ~ v u ~ 100 GeV, 
and obtain 

/i^'y ~ 2 • 10~ 4 GeV . (15) 

If we compare Eq. (|15p with the astrophysical bound in Eq. (jlip obtained from the tree level con- 
tribution to electron-Majoron coupling, assuming that the neutrino Yukawa couplings involved with 
Majoron emission and neutrino mass generation are of the same order, we obtain two limits: 

y>20keV, h^<W. (16) 



On the other hand, if we compare Eq. (|15p with Eq. (112p . the bound on the electron-Majoron coupling 
from the loop process, we obtain 

y>40eV, hW < 5 • 10 3 . (17) 

The upper bound on h^ N ' is not particularly useful if the theory has to be perturbative. However, 
it is interesting to observe that while the astrophysical bounds in Eqs. (jlip and (|12p are upper bounds 
on a combination of the neutrino Yukawa coupling and the ^-violating VEV, finally we obtain lower 
bounds on the latter — see Eqs. (]16p and (|17p — if we assume some reasonable value of the light 
neutrino mass. This is because Eqs. , (fTUj) and (HU) imply an order-of- magnitude relation between 
Majoron coupling and neutrino mass: 

tree M \ ™em» loop 1 Mg m e m v 

9eeJ ~ g 2 vlM 2 s V ' 9eeJ 16vr 2 v 2 M 2 V ' { } 

where, we recall from Eq. (|10p that M$ is either or Mz, whichever is larger. Such a proportionality 
between neutrino mass and Majoron-electron coupling occurs also in other singlet Majoron models [19J. 

It should be noted that the electron-Majoron coupling in Eqs. ([9]) and (|10p contains (h^) 2 . In 
writing this, we have suppressed the generation indices. More explicitly, the combination that actually 
appears is h>^ . The same combination is constrained from the electron-neutrino mass, whose 
expression appears in Eq. (|14p . This implies further constraints on other combinations, involving 
different charged leptons, through neutrino mixing parameters [20j . 

Now we turn our attention to the A and A' couplings in Eqs. ([3]) and and see what information 
on we get from them. Using the neutrino mass constraint in Eq. (|15p . the dimensionless prefactor 
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hWv/n~ 1(T 6 , for fj, ~ 100 GeV, provides sufficient suppression to A and A' couplings, in addition to 
those coming from charged lepton (or, down-type quarks) Yukawa couplings, to meet all experimental 
constraints |21H25| . As a result, we can keep the matrix elements to be all order unity. We 
note that, unlike the trilinear A or A' couplings, the bilinear /j, a parameters do not pick up the extra 
suppression from charged lepton Yukawa couplings, and we may expect that the corresponding bilinear 
soft terms are not suppressed either. This observation helps us to face an important question at this 
stage: how do we produce an acceptably large second neutrino mass eigenvalue? This could be induced 
by Grossman-Haber loops [26] which contribute to neutrino mass through the /^-violating bilinear soft 
terms. In these loops there are gauge couplings at the neutrino vertices, and there are two types of 
/^-violating interactions (e.g. slepton-Higgs or neutrino-neutralino mixing) inside the loop giving rise 
to AC = 2 interactions. Addition of these loops to Eqs. (|13p and (|14p breaks the rank one structure 
of the mass matrix, but one eigenvalue still remains zero. This is very much consistent with the 
neutrino oscillation data, which do not any way compel us to consider a non-vanishing third mass 
eigenvalue. The generation of the latter requires the relevant A or A' couplings to be ~ (10 -3 — 10 ) 
for superparticle masses of order 100 GeV, but in our scenario these couplings are significantly more 
suppressed (For a detailed discussion of how RPV couplings generate neutrino masses and mixing, 
see, for example, [27-33] ) . 



5 Scenario with gauged lepton number 

Eq. ([2]) can also be interpreted as the superpotential of a model where the lepton number symmetry is 
gauged. Of course, lepton number is anomalous, but the combination B — C is not. So, as a simplest 
example, we can think of a model where the gauge symmetry consists, apart from the standard model 
gauge group, of another factor of U(l)g_£. This is the same as the model presented in Ref. [34J, 
where a different combination of the weak hypercharge and B — C had been used to denote this extra 
symmetry. 

Without any loss of generality, we assume that only one singlet sneutrino acquires VEV, and to 
avoid confusion with the global symmetry case, we denote this VEV by vr. There will be no Majoron 
in this case: the Goldstone boson will be eaten up by the extra neutral gauge boson that is present in 
this model. The strength of -R-parity violation will be related to the strength of this new gauge force. 
It has been shown [33] in the context of an I?6-inspired model that 

4 25 
M 2 Z , = - tan 2 9 W M^ + — g a v\ . (19) 

Using the current experimental lower limit Mz> > 900 GeV [35] from the Tevatron pp-collider, we 
obtain the limit 

u fl >1.7TeV. (20) 

Although we cannot use the astrophysical bounds for this model since there is no Majoron, the 
neutrino mass formula given in Eq. (|14p still holds, where V should be read as vr. Using m v = 0.1 eV 
as before and putting the lower bound on vr from Eq. (I20j) , we obtain 

M x ~ (/i (A °) 2/3 x 4500 TeV. (21) 

Unlike in the previous example with global lepton number symmetry, the elements of the Yukawa 
matrix h^ N ' will now have to be very small in order to keep the neutralino mass My in the phe- 
nomenologically interesting range of a few hundred GeV to a TeV. 
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6 Conclusion 



We have done a few new things in this paper. In spite of the existence of a vast literature on 
spontaneous i?-parity violation, explicit expressions of lepton number violating couplings in terms of 
the sneutrino VEV and neutrino Yukawa couplings were somehow missing. Also, only bilinear RPV 
terms were discussed in the context of Majoron models so far. As we have shown, trilinear RPV terms 
would be generated too. We displayed them in Eqs. ([3]), flU) and ([5]). These might be particularly 
useful if one attempts to construct some flavor models relating different entries of the matrix, 
which can comfortably be of order unity. Note that using 10 parameters we can predict 39 -R-parity 
violating couplings. The other new thing that we presented is an explicit estimate of the bounds on the 
singlet sneutrino VEV and the generic fcW by using the astrophysical constraint on Majoron-electron 
coupling and the knowledge of the neutrino mass: see Eqs. (|16p and (|17p . 

We have also noted that when lepton number is gauged, the non-observation of any additional 
gauge boson in any collider experiment puts a strong bound on the singlet sneutrino VEV, which in 
turn demands that entries of the matrix have to be small to keep the neutralino masses in the 
accessible range. 
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